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1 Introduction 

Phenomena beyond the standard model (SM) of particle physics can manifest themselves di¬ 
rectly, via the production of new particles, or indirectly, by affecting the production and decay 
of SM particles. Analyses of flavor-changing neutral current (FCNC) decays are particularly 
sensitive to the effect of new physics, since such decays are highly suppressed in the SM. The 
FCNC decay, (K*° indicates the K*(892)°, and charge-conjugate states are im¬ 

plied for all particles unless stated otherwise), provides many opportunities to search for new 
phenomena. In addition to the branching fraction, other properties of the decay can be mea¬ 
sured, including the forward-backward asymmetry of the muons, A^q, and the longitudinal 
polarization fraction of the K*°, Fl. To better understand this decay, these quantities can be 
measured as a function of the dimuon invariant mass squared {(f')- New physics may modify 
any of these quantities IITI-ITtI relative to their SM values IITl lT8ti24l . While previous measure¬ 
ments by BaBar, Belle, CDF, LHCb, and CMS are consistent with the SM Il25lj29l . they are still 
statistically limited, and more precise measurements offer the possibility to uncover physics 
beyond the SM. 

In this Letter, we present measurements of ApB, Fp, and the differential branching fraction 
dB /from B'^ ^ decays, using data collected from pp collisions at the CERN 

LHC by the CMS experiment at a center-of-mass energy of 8 TeV. The data correspond to an 
integrated luminosity of 20.5 ± 0.5 fb^^ [30l. The K*° is reconstructed through its decay to 
K+ 71^, and the B° is reconstructed by fitting the two identified muon tracks and the two hadron 
tracks to a common vertex. The values of ApB and Fp are measured by fitting the distribution 
of events as a function of two angular variables: the angle between the positively charged 
muon and the B*^ in the dimuon rest frame, and the angle between the K+ and the B° in the 
K*o j- 0 st frame. All measurements are performed in bins from 1 to 19GeV^. The q^ bins 
8.68 < q^ < 10.09 GeV^ and 12.90 < q^ < 14.18 GeV^, corresponding to the B° ^ ]/xpK*^ and 
B*^ —> decays {ip' refers to the tp{2S)), respectively, are used to validate the analysis. The 

former is also used to normalize the differential branching fraction. 


2 CMS detector 

A detailed description of the GMS detector, together with a definition of the coordinate sys¬ 
tem used and the standard kinematic variables, can be found in Ref. 1311 . The main detector 
components used in this analysis are the silicon tracker and the muon detection systems. The 
silicon tracker, located in the 3.8 T field of a superconducting solenoid, consists of three pixel 
layers and ten strip layers (four of which have a stereo view) in the barrel region accompanied 
by similar endcap pixel and strip detectors on each side that extend coverage out to \ri\ < 2.5. 
For tracks with transverse momenta 1 < px < lOGeV and |?/| < 1.4, the resolutions are typi¬ 
cally 1.5% in px and 25-90 (45-150) pm in the transverse (longitudinal) impact parameter 11321 . 
Muons are measured in the range \tj\ < 2.4, with detection planes made using three technolo¬ 
gies: drift tubes, cathode strip chambers, and resistive plate chambers I133II . In addition to the 
tracker and muon detectors, GMS is equipped with electromagnetic and hadronic calorimeters 
that cover |p| <5. 

Events are selected using a two-level trigger system. The first level has specialized hardware 
processors that use information from the calorimeters and muon systems to select the most 
interesting events. A high-level trigger processor farm further decreases the event rate from 
around 90 kHz to around 400 Hz, before data storage. 
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3 Reconstruction, event selection, and efficiency 


3 Reconstruction, event selection, and efficiency 

The criteria used to select the candidate events during data taking (trigger) and after full event 
reconstruction take advantage of the fact that mesons have relatively long lifetimes and 
therefore decay on average about 1 mm from their production point. The trigger only uses 
muons to select events, while the offline selection includes the full reconstruction of all decay 
products. 

All events used in this analysis were recorded with the same trigger, requiring two identi¬ 
fied muons of opposite charge to form a vertex that is displaced from the pp collision region 
(beamspot). The beamspot position (most probable collision point) and size (the extent of the 
luminous region covering 68% of the collisions in each dimension) were continuously mea¬ 
sured through Gaussian fits to reconstructed vertices as part of the online data quality moni¬ 
toring. The trigger required each muon to have pj > 3.5GeV, Ip] < 2.2, and to pass within 2cm 
of the beam axis. The dimuon system was required to have pj > 6.9 GeV, a vertex fit proba¬ 
bility larger than 10%, and a separation of the vertex relative to the beamspot in the transverse 
plane of at least 3(7, where a includes the calculated uncertainty in the vertex position and the 
measured size of the beamspot. In addition, the cosine of the angle, in the transverse plane, be¬ 
tween the dimuon momentum vector and the vector from the beamspot to the dimuon vertex 
was required to be greater than 0.9. 

The offline reconstruction requires two muons of opposite charge and two oppositely charged 
hadrons. The muons are required to match those that triggered the event readout, and also to 
pass general muon identification requirements. These include a track matched to at least one 
muon segment (collection of hits in a muon chamber consistent with the passage of a charged 
particle), a track fit x^ per degree of freedom less than 1.8, hits in at least six tracker layers with 
at least two from the pixel detector, and a transverse (longitudinal) impact parameter with 
respect to the beamspot less than 3 cm (30 cm). The reconstructed dimuon system must also 
satisfy the same requirements that were applied in the trigger. 

The hadron tracks are required to fail the muon identification criteria, have pj > 0.8 GeV, 
and have an extrapolated distance of closest approach to the beamspot in the transverse plane 
greater than twice the sum in quadrature of the distance uncertainty and the beamspot trans¬ 
verse size. The two hadrons must have an invariant mass within 90 MeV of the accepted K*^ 
mass 13^ for either the K+tt^ or K^tt^ hypothesis. To remove contamination from (/)(1020) —> 
decays, the invariant mass of the hadron pair must be greater than 1.035 GeV when the 
charged kaon mass is assigned to both hadrons. The candidates are obtained by fitting the 
four charged tracks to a common vertex, and applying a vertex constraint to improve the res¬ 
olution of the track parameters. The B° candidates must have pj > 8 GeV, \f]\ < 2.2, vertex fit 
X^ probability larger than 10%, vertex transverse separation from the beamspot greater than 12 
times the sum in quadrature of the separation uncertainty and the beamspot transverse size, 
and cos ocxy > 0.9994, where Kxy is the angle, in the transverse plane, between the B° momen¬ 
tum vector and the line-of-flight between the beamspot and the B*^ vertex. The invariant mass m 
of the B'^ candidate must also be within 280 MeV of the accepted B*^ mass m^o IIMl for either the 
or hypothesis. The selection criteria are optimized using simulated 

signal samples (described below) and background from data using sidebands of the B'^ mass. 
After applying the selection criteria, events in which at least one candidate is found contain 
on average 1.05 candidates. A single candidate is chosen from each event based on the best B^ 
vertex fit x^- 

From the selected events, the dimuon invariant mass q and its calculated uncertainty Cq are 
used to distinguish the signal from the control samples. The control samples B*^ —> ]/tpK*^ and 
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> xp'K*^ are defined by \q — < 3crq and \q — m^/| < 3crq, respectively, where mj/^ and 

rtiq,! are the accepted masses IlMI . The average value for Cq is about 26MeV. The signal sample 
is composed of the events that are not assigned to the ]/ip and tp' samples. 

The signal sample still contains contributions from the control samples, mainly due to unre¬ 
constructed soft photons in the charmonium decay These events will have a low q value and 
fall outside the selection described above. These events will also have a low m value and there¬ 
fore they can be selectively removed using a combined selection on q and m. For q < 

{q > my^), we require \{m — m^o) — {q — my^p)\ > 160 (60)MeV. For q < m^/ {q > niq,/), we 
require |{m — ntgo) — {q — niq,i)\ > 60 (30)MeV. The requirements are set such that less than 
10% of the background events originate from the control channels. 

The four-track vertex candidate is identified as a or B^ depending on whether the K+ or 
71+ invariant mass is closest to the accepted K*® mass. The fraction of candidates assigned 
to the incorrect state is estimated from simulations to be 12-14%, depending on q^. 

The global efficiency, e, is the product of the acceptance and the combined trigger, reconstruc¬ 
tion, and selection efficiency, both of which are obtained from Monte Carlo (MC) simulations. 
The pp collisions are simulated using PYTHIA Il35]l version 6.424, the unstable particles are de¬ 
cayed by EVTGEN Il36l version 9.1 (using the default matrix element for the signal), and the 
particles are propagated through a detailed model of the detector with Geant4 Il37l . The recon¬ 
struction and selection of the generated events proceed as for data. Three simulated samples 
were created in which the B'^ was forced to decay to K**^(K+7r^)/r+^^, 

or The samples were constructed to ensure that the number and spatial 

distribution of pp collision vertices in each event match the distributions found in data. The ac¬ 
ceptance is obtained from generated events, before the particle propagation with Geant4, and 
is calculated as the fraction of events passing the single-muon requirement of Pt(f) > 3-3 GeV 
and \q{}i)\ < 2.3 relative to all events with Pt(B°) > 8GeV and |/y(B°)| < 2.2. As the ac¬ 
ceptance requirements are placed on the generated quantities, they are less restrictive than the 
final selection requirements, which are based on the reconstructed quantities, to allow for the 
effect of finite resolution. Only events passing the acceptance criteria are processed through 
the GEANT simulation, the trigger simulation, and the reconstruction software. The combined 
trigger, reconstruction, and selection efficiency is the ratio of the number of events that pass the 
trigger and selection requirements and have a reconstructed B*^ compatible with the generated 
B*^ in the event, relative to the number of events that pass the acceptance criteria. The compat¬ 
ibility of generated and reconstructed particles is enforced by requiring the reconstructed K+, 
TT^, ^+, and }i^ to have (A//)^ -|- {Acp)^ less than 0.3 (0.004) for hadrons (muons), where Aq 
and Aq) are the differences in q and f between the reconstructed and generated particles. Re¬ 
quiring all four particles in the B^ decay to be matched results in an efficiency of 99.6% (0.4% of 
the events have a correctly reconstructed B° that is not matched to a generated B°) and a purity 
of 99.5% (0.5% of the matched candidates are not a correctly reconstructed B°). Efficiencies are 
determined for both correctly tagged (the K and tt have the correct charge) and mistagged (the 
K and tt charges are reversed) candidates. 


4 Analysis method 

This analysis measures ApB, Fl/ and d6/ dq^ of the decay B° —> as a function of q^. 

Figure shows the angular observables needed to define the decay: 0 k is the angle between 
the kaon momentum and the direction opposite to the B*^ (B^) in the (K*^) rest frame, 
0/ is the angle between the positive (negative) muon momentum and the direction opposite 
to the B° (B°) in the dimuon rest frame, and (p is the angle between the plane containing the 
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4 Analysis method 


two muons and the plane containing the kaon and pion. As the extracted angular parameters 
Afb and Fi do not depend on (p and the product of the acceptance and efficiency is nearly 
constant as a function of (p, the angle (p is integrated out. Although the K+ k~ invariant mass 
must be consistent with that of a K*^, there can be a contribution from spmless (S-wave) K+ 
combinations ll24ll38] - (40] . This is parametrized with two terms: Fg, which is related to the S- 
wave fraction, and Ag, which is the interference amplitude between the S-wave and P-wave 
decays. Including this component, the angular distribution of can be written 

as 1211: 

1 d^r 

r dcos 6 k dcos 0/ dq^ 

+ 
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Figure 1: Sketch showing the definition of the angular observables 0/ (left), 0 k (middle), and (p 
(right) for the decay —)■ K*°(K+7r“)/^+f/^. 


For each q^ bin, the observables of interest are extracted from an unbinned extended maximum- 
likelihood fit to three variables: the invariant mass ni and the two angular variables 

0K and 0/. For each q^ bin, the unnormalized probability density function (PDF) has the follow¬ 
ing expression: 


PDF(m,0K,0/) =rsC 


SC(m)S“(0K,0/)eC(0K,0/) 


rM 

^I^JM S“(-0K, -0/) 6'^(0K,0/) 


( 2 ) 


where the contributions correspond to correctly tagged signal events, mistagged signal events, 
and background events. The parameters Yg and Yb are the yields of correctly tagged signal 
events and background events, respectively, and are free parameters in the fit. The parameter 
is the fraction of signal events that are mistagged and is determined from MC simula¬ 
tion. The signal mass probability functions S^{m) and S^{m) are each the sum of two Gaus¬ 
sian functions and describe the mass distribution for correctly tagged and mistagged signal 
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events, respectively. In the fit, there is one free parameter for the mass value in both sig¬ 
nal functions, while the other parameters (four Gaussian cr parameters and two fractions re¬ 
lating the contribution of each Gaussian) are obtained from MG simulation, which has been 
found to accurately reproduce the data. The function S"(0K/ Gj) describes the signal in the two- 
dimensional (2D) space of the angular observables and corresponds to Eq. Q. The combination 
B^{m) is obtained from sideband data and describes the background in the 

space of {m, 0K/ di), where the mass distribution is an exponential function and the angular dis¬ 
tributions are polynomials ranging from second to fourth degree, depending on the bin and 
the angular variable. The functions {9k, dj) and e^{6K, 9i) are the efficiencies in the 2D space 
of — 1 < COS0K < 1/ —1 < cos 9i < 1 for correctly tagged and mistagged signal events, respec¬ 
tively. The efficiency function for correctly tagged events is obtained from a fit to the 2D-binned 
efficiency from simulation and is constrained to be positive. There are 30 bins (5 in cos 0 k and 6 
in cos 0/), and the efficiency fit function is a polynomial of third degree in cos 0 k and fifth degree 
in cos 0/ (and all cross terms), for a total of 24 free parameters. This procedure does not work 
for the mistagged events because of the much smaller number of events (resulting in empty 
bins) and a more complicated efficiency. For mistagged events, the 2D efficiency is calculated 
in 5x5 bins of cos 0 k and cos 0/, and an interpolation is performed. This interpolation function 
is used to generate a new binned efficiency (in 120 x 120 bins), with all bin contents constrained 
to be normegative. The efficiency function uses this finely binned efficiency, with linear inter¬ 
polation between bins. The efficiencies for both correctly tagged and mistagged events peak at 
cos 0/ near 0 for q^ < 10 GeV^, becoming flat for larger values of q^. The efficiency for correctly 
tagged events tends to decrease with increasing cos 0 k, and for q^ > 14 GeV^ a small decrease 
is seen for cos 0 k near —1. The efficiency for mistagged events is maximal near cos 0 k = 0, with 
an increase as cos 0 k approaches -|-1 that becomes more pronounced as q^ increases. 

The fit is performed in two steps. The initial fit uses the data from the sidebands of the mass 
to obtain the B®k(0k) and B^'(0/) distributions (the signal component is absent from this fit). 
The sideband regions are 3am < \m — m-go \ < 5.5am, where am is the average mass resolution 
(?a45 MeV), obtained from fitting the MG simulation signal to a sum of two Gaussians with a 
common mean. The distributions obtained in this step are then fixed for the second step, which 
is a fit to the data over the full mass range. The free parameters in this fit are Afb, Fg, Fg, Ag, 
the parameters in B’"{m), the mass parameter in S^{m) and S^(m), and the yields Yg and Yg. 
In addition, the remaining parameters in S^(m) and S^{m) are free parameters with Gaussian 
constraints from previous fits to simulated signal events. 


The PDF in Eq. (|^ is only guaranteed to be nonnegative for particular ranges of Afb, Fg, Ag, 
and Fg. While the definition of the precise physical region is a more complicated expression, the 

approximate ranges of validity are: 0 < Fl < 1, |Afb| < 5 (1 — El)/ 0 < Fg < min 


3(l-A) 

1+3J1 


,1 


and |Ag| < Fg -|- 3 Fl (1 — Fg). In addition, the interference term Ag must vanish if either of 
the two interfering components vanish. From Ref. |24l , this constraint is implemented as 
I Ag I < ,yi2Fs(T^-FsJWR, where R is a ratio related to the S-wave and P-wave line shapes, 
estimated to be 0.89 near the mass. During the MINUIT EH minimization, penalty terms 
are introduced to ensure that parameters remain in the physical region. When assessing the 
statistical uncertainties with MlNOS BlU , the penalty terms are removed. However, a negative 
value for Eq. (j^ results in the minimizing algorithm generating a large positive jump in the 
negative log-likelihood, tending to remove the unphysical region. The results of the fit in each 
signal q^ bin are Afb, Fl, Ag, Fg, and the correctly tagged signal yield Yg. 

The differential branching fraction, d/3/dq^, is measured relative to the normalization channel 







6 


5 Systematic uncertainties 


using: 

dB (bo ^ _( y ^, \ (yN , ^n/n (b° ^ 

{l-ne^)[eC^ {l-f^)e^) Aq^ 

(3) 

where Yg and are the yields of the correctly tagged signal and normalization channels, re¬ 
spectively; Cg and are the efficiencies for the correctly tagged signal and normali z ation chan¬ 
nels, respectively; and are the mistag rates for the signal and normalization channels, re¬ 
spectively; and are the efficiencies for the mistagged signal and normalization channels, 
respectively; and B (B° ^ = 0.132% x 5.96% is the accepted branching frac¬ 

tion for the normali z ation channel 1341 , corresponding to the q^ bin Aq^ = 8.68 — 10.09 GeV^. 
The efficiencies are obtained by integrating the efficiency functions over the angular variables, 
weighted by the decay rate in Eq. Q, using the values obtained from the fit of Eq. (|^ to the 
data. 

The fit formalism and results are validated through fits to pseudo-experimental samples, MC 
simulation samples, and control channels. Additional details, including the sizes of the sys¬ 
tematic uncertainties assigned from these fits, are described in Section]^ 


5 Systematic uncertainties 

Since the efficiency is computed with simulated events, it is essential that the MC simulation 
program correctly reproduces the data, and extensive checks have been performed to verify the 
accuracy of the simulation. The systematic uncertainties associated with the efficiencies, and 
other sources of systematic uncertainty are described below and summarized in Table 

The correctness of the fit function and the procedure for measuring the variables of interest 
are verified in three ways. Eirst, a high-statistics MC sample (approximately 400 times that 
of the data) is used to verify that the fitting procedure produces results consistent with the 
input values to the simulation. This MC sample includes the full simulation of signal and con¬ 
trol channel events plus background events obtained from the PDE in Eq. Q. The discrepancy 
between the input and output values in this check is assigned as a simulation mismodeling sys¬ 
tematic uncertainty. It was also verified that fitting a sample with only mistagged events gives 
the correct results. Second, 1000 pseudo-experiments, each with the same number of events as 
the data sample, are generated in each q'^ bin using the PDE in Eq. (j^, with parameters obtained 
from the fit to the data. These are used to estimate the fit bias. Much of the observed bias is 
a consequence of the fitted parameters lying close to the boundaries of the physical region. In 
addition, the distributions of results are used to check the returned statistical uncertainty from 
the fit and are found to be consistent. Third, the high-statistics MC signal sample is divided 
into 400 subsamples and combined with background events to mimic 400 independent data 
sets of similar size to the data. Eits to these 400 samples do not reveal any additional systematic 
uncertainty. 

Because the efficiency functions are estimated from a finite number of simulated events, there 
is a corresponding statistical uncertainty in the efficiency. The efficiency functions are obtained 
from fits to simulated data. Alternatives to the default efficiency function are generated by ran¬ 
domly varying the fitted parameters within their uncertainties (including all correlations). The 
effect of these different efficiency functions on the final result is used to estimate the systematic 
uncertainty. 

The main check of the correctness of the efficiency is obtained by comparing the efficiency- 
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corrected results for the control channels with the corresponding world-average values. The 
efficiency as a function of the angular variables is checked by comparing the Fl and Afb mea¬ 
surements from the > ]/tpK*^ sample, composed of 165 000 signal events. The value of 
Fl obtained in this analysis is 0.537 ± 0.002 (stat), compared with the world-average value of 
0.571 ± 0.007 (stat+syst) Il34ll , indicating a discrepancy of 0.034, which is taken as the system¬ 
atic uncertainty for the signal measurements of Fl. For ApB/ the measured value is 0.008 ± 
0.003 (stat), compared to a SM expectation of «0. Adding an S-wave contribution in the fit 
changes the measured value of ApB by less than 0.001. From this, we conclude that the S-wave 
effects are minimal, and assign a systematic uncertainty of 0.005 for Afb- To validate that the 
simulation accurately reproduces the efficiency as a function of we measure the branch¬ 
ing ratio between two different q^ bins, namely the two control channels. The branching ratio 
result, B (B° —> IB (B° —> J/t/;K*°) = 0.479 ± 0.005, is in excellent agreement with the 

most precise reported measurement: 0.476 ± 0.014 (stat) ± 0.010 (syst) Il42l . 

The PDF used in the analysis accommodates cases in which the kaon and pion charges are 
correctly and incorrectly assigned. Both of these contributions are treated as signal. The mistag 
fraction is fixed to the value obtained from MC simulation. In the high-statistics control channel 
B*^ —)■ ]/xpK*^, the mistag fraction is allowed to float in the fit and a value of = (14.5 ± 
0.5)% is found, to be compared to the simulated value of (13.7 ± 0.1)%. The effect of this 5.8% 
difference in the mistag fraction on the measured values is taken as a systematic uncertainty. 

The systematic uncertainty associated with the functions used to model the angular distribu¬ 
tion of the background is obtained from the sum in quadrature of two uncertainties. The first 
uncertainty is evaluated by fitting the background with polynomials of one degree greater than 
used in the default analysis and taking the difference in the observables of interest between 
these two fits as the systematic uncertainty. The second uncertainty is owing to the statistical 
uncertainty in the background shape, as these shapes are fixed in the final fit. This uncertainty 
is obtained by taking the difference in quadrature between the returned statistical uncertainties 
on the parameters of interest when the background shapes are fixed and allowed to vary. In 
q^ bins where the unconstrained fit does not converge, the associated uncertainty is obtained 
from extrapolation of nearby bins. 

The mass distributions for the correctly tagged and mistagged events are each described by 
the sum of two Gaussian functions, with a common mean for all four Gaussian functions. The 
mean value is obtained from the fit to the data, while the other parameters (four cr and two 
ratios) are obtained from fits to MG-simulated events, with the uncertainty from those fits used 
as Gaussian constraints in the fits to the data. For the high-statistics control channels, it is 
possible to fit the data, while allowing some of the parameters to vary. The maximum changes 
in the measured values in the two control channel q^ bins when the parameters are varied are 
taken as the systematic uncertainty for all q^ bins. 

The q^ bins just below and above the ]/xp region may be contaminated with B*^ —> ]/ipK*^ feed¬ 
through events that are not removed by the selection criteria. A special fit in these two bins 
is made, in which an additional background term is added to the PDF. This background dis¬ 
tribution is obtained from the MG simulation and the background yield is a free parameter. 
The resulting changes in the fit parameters are used as estimates of the systematic uncertainty 
associated with this contribution. 

The effects from angular resolution in the reconstructed values for the angular variables 0 k and 
0/ are estimated by performing two fits on the same MG-simulated events. One fit uses the true 
values of the angular variables and the other fit their reconstructed values. The difference in 
the fitted parameters between the two fits is taken as an estimate of the systematic uncertainty. 
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6 Results 


The differential branching fraction has an additional systematic uncertainty of 4.6% coming 
from the uncertainty in the branching fraction of the normalization mode —)■ ]/ipK*^. 

The systematic uncertainties are measured and applied in each bin, with the total systematic 
uncertainty obtained by adding the individual contributions in quadrature. 

Table 1: Systematic uncertainty contributions for the measurements of Fl, A^b, and the branch¬ 
ing fraction for the decay —>• The values for Fl and ApB are absolute, while the 

values for the branching fraction are relative. The total uncertainty in each q^ bin is obtained 
by adding each contribution in quadrature. For each item, the range indicates the variation of 
the uncertainty in the signal q^ bins. 


Systematic uncertainty 

Fl(10-3) 

Afb(10 ^) 

dB/dq^ (%) 

Simulation mismodeling 

1-17 

0-37 

1.0-5.5 

Fit bias 

0-34 

2-42 

— 

MG statistical uncertainty 

3-10 

5-18 

0.5-2.0 

Efficiency 

34 

5 

— 

Ktt mistagging 

1-4 

0-7 

0.1-4.1 

Background distribution 

20-36 

12-31 

0.0-1.2 

Mass distribution 

3 

1 

3.2 

Feed-through background 

0-27 

0-5 

0.0-4.0 

Angular resolution 

6-24 

0-5 

0.2-2.1 

Normalization to B'^ ^ ]/tpK*^ 

— 

— 

4.6 

Total systematic uncertainty 

41-65 

18-74 

6.4-8.6 


6 Results 

The signal data, corresponding to 1430 signal events, are fit in seven disjoint q^ bins from 1 
to 19GeV^. Results are also obtained for a wide, low-q^ bin {1 < q^ < 6GeV^), where the 
theoretical uncertainties are best understood. The invariant mass distributions for 

all of the q^ signal bins, as well as the fit projections, are shown in Fig. Figure plots the 
projections of the fit and the data on the cos 0 k (top) and cos 0/ (bottom) axes for the combined 
\ow-q^ bin (left, 1 < q^ < 6GeV^) and the highest q^ bin (right, 16 < q^ < 19GeV^). The fitted 
values of signal yield, Fl, Alb, and diB/ dq^, along with their associated uncertainties, are given 
for each of the disjoint q^ regions in Table These results are also shown in Fig. along with 
two SM predictions. The fitted values for Fs are all less than 0.03, while the values for As vary 
from —0.3 to +0.3. 

The SM predictions, derived from Refs. IITSlI^ . combine two calculational techniques. In the 
\ow-q^ region, a quantum chromodynamic factorization approach Il4^ is used, which is appli¬ 
cable for q^ < Am^, where rric is the charm quark mass. In the high-q^ region, an operator prod¬ 
uct expansion in the inverse b quark mass and 1 / IImI |4^ is combined with heavy-quark 
form-factor relations |46J. This is valid above the open-charm threshold (q^ > 13.9GeV^). The 
two SM predictions shown in Fig. [^differ in the calculation of the form factors. The light-cone 
sum rules (LGSR) calculation is made at low q^ 147 1 and is extrapolated to high q^ i^|. The 
lattice gauge (Lattice) calculation of the form factors is from Ref. Il49l . Gontrolled theoretical 
predictions are not available near the ]/tp and xp' resonances. The SM predictions are in good 
agreement with the GMS experimental results, indicating no strong contribution from physics 
beyond the standard model. 
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Figure 2: The invariant mass distributions for the seven signal bins and the 

combined 1 < q^ < 6 GeV^ bin. Overlaid on each is the projection of the results for the total fit, 
as well as the three components: correctly tagged signal, mistagged signal, and background. 
The vertical bars give the statistical uncertainties, the horizontal bars the bin widths. 
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7 Summary 




Figure 3: Data and fit results for 1 < < 6 GeV^ (left) and 16 < q^ < 19 GeV^ (right), projected 

onto the cos 0 k (top), and cos 0/ axis (bottom). The fit results show the total fit, as well as 
the three components: correctly tagged signal, mistagged signal, and background. The vertical 
bars give the statistical uncertainties, the horizontal bars the bin widths. 


The results described are combined with previous GMS measurements, obtained from an in¬ 
dependent data sample collected at ^/s = 7TeV Il29ll . The systematic uncertainties associated 
with the efficiency, Krc mistagging, mass distribution, angular resolution, and the —)■ 

branching fraction are assumed to be fully correlated between the two samples, with the re¬ 
maining uncertainties assumed to be uncorrelated. To combine the results from the 7TeV and 
8 TeV data, the uncorrelated systematic uncertainties are combined in quadrature with the sta¬ 
tistical uncertainties. To account for the asymmetric uncertainties, the linear variance method 
from Ref. ||50]| is used to average the 7 TeV and 8 TeV measurements, as well as to average the 
two q^ bins covering 4.30 to 8.68 GeV^, which was a single bin in the 7 TeV analysis. After the 
combination, the correlated systematic uncertainties are added in quadrature. The combined 
GMS measurements of Afb, Fl/ arid the differential branching fraction versus q^ are compared 
to previous measurements ll26ti29llMll52]| in Fig.|^ The GMS measurements are consistent with 
the other results, with comparable or higher precision. Table [^provides a comparison of the 
measured quantities in the low dimuon invariant mass region: 1 < < 6 GeV^, as well as the 

corresponding theoretical calculations. 


7 Summary 

Using pp collision data recorded at = 8 TeV with the GMS detector at the LHG, correspond¬ 
ing to an integrated luminosity of 20.5 fb^^, an angular analysis has been carried out on the 
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Figure 4: Measured values of Fi, Afb/ and dB/dcj^ versus for The statistical 

uncertainty is shown by the inner vertical bars, while the outer vertical bars give the total 
uncertainty. The horizontal bars show the bin widths. The vertical shaded regions correspond 
to the J/tp and ip' resonances. The other shaded regions show the two SM predictions after rate 
averaging across the bins to provide a direct comparison to the data. Controlled theoretical 
predictions are not available near the ]/xp and xp' resonances. 
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7 Summary 




Figure 5: Measured values of Fl, Afb, and <\BI versus for from CMS 

(combination of the 7TeV Il29l results and this analysis). Belle l26l . CDF W7\ 1^ . BaBar Il52ll . 
and LHCb |28|. The CMS and LHCb results are from B^ —>• decays. The remaining 

experiments add the corresponding B+ decay, and the BaBar and Belle experiments also include 
the dielectron mode. The vertical bars give the total uncertainty. The horizontal bars show the 
bin widths. The horizontal positions of the data points are staggered to improve legibility. The 
vertical shaded regions correspond to the J/t/? and xp' resonances. 
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Table 2: The measured values of signal yield (including both correctly tagged and mistagged 
events), Fl, Aps, and differential branching fraction for the decay —>• in bins of 

The first uncertainty is statistical and the second (when present) is systematic. The bin ranges 
are selected to allow comparisons to previous measurements. 



Signal 

Fl 

ApB 

dB/dq^ 

(GeV^) 

yield 



(10-®GeV-2) 

1.00-2.00 

84 ±11 

0.64+ °;J°± 0.07 

-0.27 + °;“ ±0.07 

4.6 ±0.7 ±0.3 

2.00-4.30 

145 ± 16 

0.80 ± 0.08 ± 0.06 

-0.12l°-“±0.05 

3.3 ±0.5 ±0.2 

4.30-6.00 

117 ± 15 

0.62 + ° “ ±0.07 

0.01 ±0.15 ±0.03 

3.4 ± 0.5 ± 0.3 

6.00-8.68 

254 ± 21 

0.50 ± 0.06 ± 0.06 

0.03 ±0.10 ±0.02 

4.7 ±0.4 ±0.3 

10.09-12.86 

362 ± 25 

0.39 ± 0.05 ± 0.04 

0.16 ±0.06 ±0.01 

6.2 ± 0.4 ± 0.5 

14.18-16.00 

225 ± 18 

0.48 0.04 

0.39 0.01 

6.7 ±0.6 ±0.5 

16.00-19.00 

239 ± 18 

0.38 0.04 

0.35 ± 0.07 ± 0.01 

4.2 ± 0.3 ± 0.3 


Table 3: Measurements from CMS (the 7 TeV results (291/ this work for 8 TeV, and the com¬ 
bination), LHCb ESi . BaBar (52l . CDF HZlISTl, and Belle (2611 of fp/ ApB, and dB/dq^ in the 
region 1 < q^ < 6GeV^ for the decay B*^ —> The CMS and LHCb results are from 

B*^ —)■ decays. The remaining experiments add the corresponding B^ decay, and the 

BaBar and Belle experiments also include the dielectron mode. The first uncertainty is statisti¬ 
cal and the second is systematic. For the combined CMS results, only the total uncertainty is 
reported. The two SM predictions are also given. 


Experiment 

Fl 

ApB 

dB/dq^ (10-®GeV-2) 

GMS (7 TeV) 

0.68 ±0.10 ±0.02 

-0.07 ±0.12 ±0.01 

4.4 ± 0.6 ± 0.4 

GMS (8 TeV, this analysis) 

0.73 ± 0.05 ± 0.04 

-0.16;° “ ±0.05 

3.6 ±0.3 ±0.2 

GMS (7 TeV+ 8 TeV) 

0.72 ± 0.06 

-0.12 ±0.08 

3.8 ± 0.4 

LHGb 

0.65 lg;g^± 0.03 

-0.17 ±0.06 ±0.01 

3.4±0.3;g;| 

BaBar 

— 

— 

4.i;+J±0.1 

GDF 

Belle 

0.69 + °;“ ±0.08 
0.67 ± 0.23 ± 0.05 

0 . 29 ;°;“ ±0.07 
0.26;°;“ ±0.07 

3.2 ±1.1 ±0.3 
3.0;°|±0.2 

SM (LGSR) 

0 79 + ° 

—0 02 + ° °^ 

- 0.02 

4.6;;^ 

SM (Lattice) 

n 70 + 0.08 
^•^‘^-O.IO 

—0 03 + ° °'^ 

- 0.03 

3-8 ;i:^ 


decay B*^ —>• The data used for this analysis include 1430 signal decays. For each 

bin of the dimuon invariant mass squared {q^), unbinned maximum-likelihood fits were per¬ 
formed to the distributions of the invariant mass and two decay angles, to obtain 

values of the forward-backward asymmetry of the muons, ApB, the fraction of longitudinal po¬ 
larization of the K*^, fp/ and the differential branching fraction, dB/dq^. The results are among 
the most precise to date and are consistent with standard model predictions and previous mea¬ 
surements. 
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